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Preface 


As  part  of  this  project,  an  apparatus  was  designed  to 
study  the  free  convection  across  an  air  layer,  enclosed  in 
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tance  in  conducting  the  experiment  and  Mr.  Allen  Hause  of 
Aero-Propulsion  Laboratory  for  his  help  in  the  use  of  a 
densitometer.  I  extend  my  thanks  and  deepest  admiration 
to  Dr.  James  E.  Hitchcock  of  the  Department  of  Aeronautics 
and  Astronautics,  Air  Force  Institute  of  Technology.  His 
guidance  and  professionalism,  both  as  a  professor  and  an 
advisor,  were  greatly  appreciated.  Finally,  I  wish  to 
acknowledge  my  gratefulness  to  my  wife  and  son  for  their 
cheerful  patience  during  this  project. 

Mohammad  Khalid  Khan 
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Abstract 


A  small  experimental  apparatus  was  designed  for  the 
measurement  of  heat  flux  across  an  enclosed,  rectangular 
air  cavity,  which  could  be  placed  at  any  angle  of  inclina¬ 
tion.  The  heat  flux  evaluation  was  done  with  the  help  of 
an  interferometer.  Heat  flux  was  determined  for  five 
aspect  ratios,  each  placed  at  four  different  angles  and 
two  different  temperature  differences  for  each  setting. 

The  average  Nusselt  number  was  found  to  be  primarily  a 
function  of  Grashof  number  or  Rayleigh  number.  Three  cor¬ 
relation  equations  developed  by  earlier  studies  are  fitted 
to  the  data.  The  correlation  equations  are  found  to 
approximate  the  heat  flux  within  +  11%,  up  to  a  Grashof 
number  of  1.2  million.  - 


I.  INTRODUCTION 


Free  convection  between  two  vertical  or  inclined 
plates  has  recently  become  very  important  to  the  designers 
of  solar  collectors.  Smaller  collector  areas  can  be  used 
if  heat  loss  can  be  reduced  from  the  absorber  plate  through 
the  cover  plates.  In  other  applications,  the  convection 
arising  in  such  layers  is  utilized  to  transport  heat.  The 
problem  of  predicting  the  amount  of  heat  transported 
through  a  vertical  air  layer  has  found  considerable  atten¬ 
tion  in  the  past.  One  of  the  first  publications  on  the 
subject  was  written  by  W.  Nusselt  in  1909  (Reference  1). 

It  was  followed  by  a  series  of  other  papers,  the  most 
extensive  ones  by  Mull  and  Reiher  (Reference  2)  and  by 
DeGraaf  and  Van  der  Held  (Reference  3)-  The  primary  aim 
of  these  papers  was  to  measure  the  total  amount  of  heat 
which  is  transported  from  the  hot  plate  to  the  cold  plate 
in  a  vertical  geometry.  The  results  were  presented  in  the 
dimensionless  form: 


Nu  =  C(Gr)m  (£)n 

where  the  mean  Nusselt  number  Nu,  and  the  Grashof  number 
Gr,  are  based  on  the  width  of  the  air  layer  L,  and  H  is 
the  height  of  the  air  layer.  The  exponents  m  and  n  acquire 
different  values  for  laminar  and  turbulent  flow. 

The  problem  of  predicting  the  amount  of  heat  trans¬ 
fer  between  the  hot  and  the  cold  plate  of  a  horizontal, 


rectangular  geometry  has  also  been  studied  to  some  extent 
(Reference  4) .  Important  investigations  of  free  convection 
across  inclined  air  layers,  with  various  angles  of  incli¬ 
nation,  were  carried  out  by  DeGraaf  and  Van  der  Held 
(Reference  5)  and  by  K.  G.  T.  Hollands,  T.  E.  (Jnny, 

G.  D.  Raithby  and  L.  Konicek  (References  6  and  7) • 

However  the  authors  have  limited  their  studies  to  Rayleigh 
numbers  ranging  from  1000  to  100,000.  The  more  recent  work 
by  K.  G.  T.  Hollands  and  others  (Reference  6)  has  provided 
a  correlation  of  the  form 

Nu  =  Nu(Ra  cos  <j> ) 

where  Nu  is  given  as  a  function  of  Rayleigh  number,  Ra  = 

Gr  •  Pr,  and  the  angle  of  inclination,  <f> . 

The  present  study  has  been  carried  out  with  the 
objective  of  confirming  the  accuracy  of  the  correlation 
equations  developed  for  the  average  Nusselt  number.  In 
particular,  it  was  desired  to  check  their  applicability 
for  Grashof  numbers  greater  than  100,000.  For  free  con¬ 
vective  heat  transfer  across  inclined  air  layers,  the 
correlation  equation  developed  by  K.  G.  T.  Hollands  and 
others  (Reference  6)  has  been  used.  For  heat  transfer 
across  vertical  air  layers,  the  correlation  equations 
developed  by  E.  R.  G.  Eckert  and  Walter  0.  Carlson  (Ref¬ 
erence  8)  have  been  used. 

In  this  research  free  convection  has  been  studied 
for  five  different  aspect  ratios,  H/L  =  11.295,  7-867, 


3  -  933  s  1.967,  0.929,  and  at  four  different  angles,  4>  =  0, 

30,  60,  and  90  degrees.  In  any  of  the  inclined  positions, 
the  hot  surface  remained  below.  The  study  has  been  carried 
out  for  two  different  temperature  differences  for  each 
geometry  and  angle.  Gas  temperature  measurements  were  made 
optically  with  the  help  of  a  Mach-Zehnder  Interferometer. 

An  apparatus  was  designed  for  the  purpose  of  this  experi¬ 
ment  which  consists  of  two  copper  plates,  each  one  foot 
long,  six  inches  wide,  and  a  quarter  inch  thick.  The 
hollow  plates  were  kept  at  a  constant  temperature  with  the 
help  of  running  cold  or  hot  water  through  the  respective 
plates.  Wood  blocks  separate  the  copper  plates,  and  opti¬ 
cal  glass  plates  seal  the  ends  of  the  rectangular  enclosure. 
A  detailed  description  of  the  apparatus  is  given  in  Section 
II.  The  experimental  procedures  and  the  data  reduction 
techniques  are  included  in  Section  III. 

It  was  proved  that  Eckert's  correlation  for  free 
convective  heat  transfer  between  two  vertical  plates  and 
Hollands'  correlation  for  free  convection  between  two 
inclined  or  horizontal  plates  hold  good  at  Rayleigh  numbers 
as  high  as  one  milli  m.  A  general  pattern  for  the  distri¬ 
bution  of  the  local  Nusselt  numbers  was  observed;  but  since 
there  was  a  lot  of  scatter  in  the  data,  these  results  are 
not  as  dependable  as  the  mean  Nusselt  number  data.  A 
detailed  discussion  about  the  results  is  included  in 


Section  IV. 


II.  APPARATUS 


Design  Concept  and  Constraints 

The  apparatus  was  designed  to  be  compatible  with  the 
23  inches  wide  test  section  of  a  Mach-Zehnder  interfero¬ 
meter.  Also,  because  the  interferometer  has  7-inch  dia¬ 
meter  mirrors,  the  maximum  length  of  the  rectangular  cavity 
was  constrained  to  6  inches.  The  layout  of  the  apparatus 
is  given  in  Figure  1(a).  The  apparatus  can  be  seen  in 
detail  in  Figure  1(b)  and  Figure  2.  The  apparatus  was 
designed  so  that  the  cavity  could  be  rotated  to  any  desired 
angle  from  the  horizontal.  Both  ends  of  the  cavity  have 
been  manufactured  from  hard  wood,  so  that  adiabatic  end¬ 
boundary  conditions  could  be  approached.  Interchangeable 
blocks  of  various  size  determine  the  aspect  ratio  of  the 
cavity.  The  optical  path  length,  12  inches,  was  chosen  to 
produce  reasonable  interference  fringe  deflections  for 
temperature  differences  up  to  35°F  between  the  isothermal 
plates.  Optical  glass  plates  at  each  end  complete  the 
enclosure.  The  temperatures  of  the  hot  and  the  cold 
copper  plates  are  maintained  constant  with  the  help  of 
running  constant  temperature  hot  or  cold  water  through 
them,  respectively. 
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An  unsealed  three-dimensional  drawing  of  the  apparatus 


is  shown  in  Figure  2. 

Copper  plates.  Two  copper  plates,  each  one  foot  long, 
six  inches  wide  and  one-quarter  inch  thick,  have  been  used 
for  the  hot  and  the  cold  surfaces  of  the  rectangular  cavity 
to  be  studied.  The  copper  plates  have  channels  built  into 
each  plate  to  form  a  hollow  body.  Each  plate  has  a  flow 
channel  entrance  and  exit  at  diagonally  opposite  corners. 
These  have  been  formed  by  drilling  holes  through  the  outer 
skin  of  the  plates  and  then  welding  a  5/8"  external  dia¬ 
meter,  1/2"  long  copper  pipe  to  the  plates  over  the  holes. 

A  sketch  of  the  plate  is  presented  in  Figure  3.  To  mini¬ 
mize  external  heat  loss,  the  hot  plate  has  been  covered 
externally  with  loose  glass-wool. 


Figure  3*  Top  View  of  the  Copper  Plate. 
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Thermocouples .  Copper-constantan  thermocouples  were 
soldered  at  three  points  into  drilled  cavities  on  the 
external  surface  of  each  plate,  as  shown  in  Figure  3*  A 
negligible  temperature  difference  between  the  inner  and 
outer  surfaces  of  the  plates  can  be  assumed  because  of 
the  high  thermal  conductivity  of  copper. 

Wooden  blocks.  The  one-foot  long  sides  of  the  copper 
plates  were  inserted  into  two  grooved  wooden  blocks  at  both 
the  sides,  thus  forming  a  rectangular  cavity  with  copper 
plates  on  two  sides  and  wooden  blocks  on  the  other  two 
sides.  Figure  4(a)  shows  an  assembled  cross-sectional 
view  of  the  5*9"  x  1.5"  geometry.  The  other  three  geome¬ 
tries  having  a  5.9"  height  look  similar.  Figure  4(b)  shows 
the  cross-sectional  view  of  the  2.786"  x  3.00"  geometry. 


Figure  4 

a)  End  View  of  Rectangular  Cavity,  5-9"  x  1.5". 

b)  End  View  of  Rectangular  Cavity,  2.786"  x  3.00". 
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ft 

The  five  different  cross-sections  thus  formed  are  listed 
below. 


H 

L 

Aspect  Ratio 

5-93" 

0.525" 

11.295 

5.9" 

0.75" 

7.867 

5.9" 

1.5" 

3-933 

5.9" 

3. 00" 

1.967 

2.786" 

3. 00" 

0.929 

The  grooves  where  the  copper  plates  fit  into  the  wooden 
blocks  have  a  lining  of  soft,  spongy  material  so  that  air 
leaks  from  the  rectangular  cavity  are  minimal.  The  two 
wooden  blocks  are  held  together  by  four  wooden  strips  on 
the  outer  side  of  the  copper  plates  so  that  the  system  is 
rigid. 

Glass  plates.  Both  open  ends  of  the  rectangular 
cavity  formed  by  the  copper  plates  and  the  wooden  blocks 
are  closed  with  the  help  of  1/4"  thick,  interferometer 
quality  glass  plates.  These  glass  plates  are  of  high 
quality,  no  air  bubbles  or  other  constructional  deforma¬ 
tions,  and  have  uniform  thickness.  The  glass  plates  are 
held  in  place  by  four  metallic  clamps  which  are  screwed 
into  the  wooden  blocks.  See  Figure  1(b).  Soft  material 
lining  is  placed  between  the  glass  plates  and  the  wooden 
blocks  to  minimize  air  leaks.  Care  was  taken  while  plac¬ 
ing  the  lining  so  that  the  view  of  the  cavity  was  not 
obstructed. 

Circular  end  supports.  The  assembly  of  copper  plates 
wooden  blocks,  and  glass  plates  was  supported  at  both  the 


front  and  the  rear  ends  by  two  circular  supports,  1.5" 
thick.  As  shown  in  Figure  5,  these  are  14.5"  diameter 
with  a  ?.5"  x  10"  rectangle  cut  out  from  their  centers. 

One  of  the  circular  supports  is  graduated  into  360°  and 
marked  at  15°  increments .  The  0°  and  180°  graduations  are 
exactly  perpendicular  to  the  copper  plates,  and  the  90°  and 
270°  graduations  are  exactly  parallel  to  the  copper  plates. 


Figure  5.  Front  View  of  the  Circular  End-Support. 

Base  support.  The  base-support,  as  shown  in  Figure  1, 
has  two  semi-circular  mounts.  The  circular  end-supports 
rest  on  top  of  these  semi-circular  mounts.  Thus  by 


rotating  the  circular  end-supports  in  the  semi-circular 
base  support,  the  copper  plates  can  be  placed  at  any 
desired  angle.  An  arrow  mark  has  been  drawn  on  the  semi¬ 
circular  base-support  in  such  a  way  that  when  0°  graduation 
of  the  circular  end-support  is  in  line  with  this  mark,  the 
copper  plates  are  exactly  horizontal.  The  circular  end- 
supports  can  be  rotated  to  any  other  angle  by  simply 
aligning  that  particular  angle  on  the  circular  end-support 
with  the  mark  on  the  base-support.  The  whole  apparatus 
once  assembled  was  placed  on  top  of  an  adjustable  tripod, 
as  shown  in  Figure  1(a),  and  can  be  rotated  and  moved  up 
and  down.  These  adjustments  were  made  to  align  the  axis 
of  the  rectangular  cavity  parallel  to  and  centered  in  the 
interferometer  light  beam. 

Other  Accessories 

Temperature  measurement  apparatus .  The  thermocouples 
were  connected  to  a  Rubicon  Potentiometer,  through  a  multi¬ 
terminal  junction  box.  The  operation  of  the  potentiometer 
is  explained  in  Reference  12.  The  cold  junction  of  the 
thermocouples  was  kept  at  32°F  with  the  help  of  an  ice- 
bath.  The  millivolts  reading  of  the  potentiometer  can  be 
converted  to  temperature  reading  through  standard  copper- 
constantan  thermocouple  conversion  table. 

Precision  constant  temperature  circulating  system. 
Constant  temperature  hot  water  was  provided  by  a  constant 
temperature  circulating  system.  The  system  consists  of  a 
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2  1/2  gallon  tank,  electric  heating  coils,  a  stirrer  and  an 
electric  pump.  The  device  also  has  an  "Adjustable  Micro- 
Set  Thermoregulator",  which  can  keep  the  temperature  of  the 
circulating  hot  water  constant  within  +  0.02°P.  The  maxi¬ 
mum  rated  temperature  is  210°F.  The  temperature  of  the 
water  in  the  bath  can  be  noted  by  inserting  a  thermometer 
into  the  bath  through  a  hole  meant  for  that  purpose.  Com¬ 
plete  operating  instructions  are  available  in  Reference 
A  picture  of  the  system  is  provided,  Figure  6. 

Cold  tap-water  was  circulated  through  the  cold  plate. 


Figure  6.  Precision  Constant  Temperature 
Circulating  System. 
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Interferometer .  A  Mach-Zehnder  interferometer  was 


used  to  measure  the  gas  temperatures  in  the  rectangular 
cavity.  A  schematic  diagram  of  this  instrument  is  given 
in  Figure  7.  Monochromatic  light  from  source  S  is  made 
parallel  by  a  parabolic  mirror  M  and  is  divided  by  a  beam 
splitter  S^,  which  is  a  semi-transparent  mirror.  The 
reflected  beam  proceeds  to  plane  mirror  Mg.  The  two  beams 
are  combined  by  splitter  plate  Sg,  which  is  identical  with 
.  The  test  section  is  between  Mg  and  Sg.  When  the 
instrument  is  properly  adjusted,  these  beams  produce  inter 
ference  fringes  which  are  then  focused  on  the  screen  or 
camera  I  by  parabolic  mirror  M^ .  A  detailed  discussion 
about  this  instrument  is  available  in  References  10  and  11 
A  polaroid  land  camera  was  used  to  photograph  the  inter¬ 
ference  fringes . 


Figure  7.  A  schematic  diagram  of  the  Mach  Zehnder  interferometer. 


Compensating  glass.  Since  the  light  beam  between  and 
S,  has  to  pass  through  two  1/4"  thick  glass  plates  of  the 
test  apparatus,  it  is  therefore  compensated  by  placing  a 
1/2"  thick  glass  plate  in  the  reference  beam  between 
and  S^.  This  is  also  a  high  quality  glass  plate  with  no 
air  bubbles. 

Microcomparator .  This  device  was  used  to  measure  the 
distances  on  an  interf erogram  (picture  of  the  fringes  from 
an  interferometer).  As  shown  in  Figure  8,  it  is  a  micro¬ 
scope  mounted  on  a  traversing  mechanism.  A  distance  on  an 
interf erogram  can  be  measured  to  ten  thousandths  of  an  inch 
with  this  device. 
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III.  EXPERIMENTAL  PROCEDURE  AND 
DATA  REDUCTION  TECHNIQUE 

Experimental  Procedure 

The  apparatus,  once  placed  in  the  test  section,  was 
adjusted  on  the  tripod  to  align  the  axis  of  the  rectangular 
cavity  with  the  interferometer  light  beam.  The  copper 
plates  were  placed  vertically  for  the  first  test  at  each 
aspect  ratio,  because  that  was  the  only  position  when  the 
complete  cross-section  could  be  seen  on  the  screen.  The 
interferometer  was  adjusted  to  give  a  set  of  finite  fringes 
perpendicular  to  the  copper  plates,  and  a  picture  of 
fringes  was  taken  in  this  position  with  a  polaroid  camera. 
Cold  tap  water  was  then  circulated  through  the  cold  plate. 
The  electrical  resistance  heaters  in  the  constant  tempera¬ 
ture  bath  were  turned  on,  and  heated  water  was  circulated 
through  the  hot  plate.  The  desired  steady-state  thermal 
condition  was  obtained  by  adjusting  the  thermoregulator 
on  the  constant  temperature  bath.  A  picture  of  the  result¬ 
ing  deflected  finite  fringe  pattern  was  then  taken.  Output 
of  the  six  copper-constantan  thermocouples  was  recorded  for 
the  accurate  determination  of  the  copper  plate  temperatures 
The  thermoregulator  was  next  set  to  a  higher  temperature, 
and  the  forementioned  procedure  was  repeated.  The  appara¬ 
tus  was  allowed  to  cool  to  room  temperature  before  it  was 
rotated  to  another  angle  of  inclination.  The  interfero¬ 
meter  was  again  adjusted  to  give  a  set  of  finite  fringes 
perpendicular  to  the  copper  plates  and  the  forementioned 
procedure  was  repeated. 
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Data  Reduction  Technique 

Finite  Fringe  Method.  The  interferometer  was  adjusted 
so  that  alternate  dark  and  light  fringes  were  obtained  as 
shown  in  Figure  9(a).  This  was  done  by  rotating  the  split¬ 
ter  lens  so  that  the  two  coherent  wave  fronts  were  combined 
at  a  small  angle.  The  light  was  then  in  phase  at  the  line 
of  intersection  and  increasingly  out  of  phase  away  from  the 
line  of  intersection.  Bright  fringes  occurred  where  the 
light  was  out  of  phase  by  0,  1,  2,  3>  .  ,  wavelengths, 

and  dark  fringes  occurred  at  the  points  where  the  light  was 
out  of  phase  by  1/2,  1  1/2,  2  1/2,  ...  ,  wavelengths.  Once 
the  temperatures  of  the  copper  plates  were  varied,  density 
variations  occurred  in  the  air  layer  enclosed  between  the 
plates,  and  the  fringes  in  the  test  section  were  displaced 
as  shown  in  Figure  9(b).  It  will  be  proved  in  Appendix  A 
that  if  T^(°R)  is  the  temperature  of  the  hot  plate  then 


T  =  - - - 

X2  1  +  e1C1T1 


where  Tg  =  temperature  at  a  distance  y^  from  the  plate 
(°R) 

y^  =  distance  of  the  first  fringe  from  the  wall  (ft) 


=  displacement  of  the  first  fringe  from  its 
initial  position  (inches) 
d  =  distance  between  fringes  (inches) 
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(Ao)(R) 

C1  =  (L}(Kdg)(P) 


Xq  =  wavelength  of  light 


K 


DG 


Dale-Gladstone  constant 


x 

o 


L 

P 


Ll  lbm 

4.91  X  10"4 

length  of  test  section  (one  foot) 
atmospheric  pressure  (lb^/ft2) 


f t-lbf 

R  =  gas  constant  for  air  =  93.35  — — - 

lb  °R 
m 


Hot  Cold 
Plate  Plate 


(a) 


(b) 


Figure  9 

a)  Initial  fringe  pattern  when  both  plates  are  at 
equal  temperature. 

b)  Fringe  pattern  when  fringes  are  displaced  because 
of  temperature  difference  between  the  plates. 
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In  a  similar  way 


T3  ~  1  +  e2CiTi 

where  =  temperature  at  a  distance  y^  from  the  plate  (°R) 

y2  =  distance  of  the  second  fringe  from  the  plate 
(ft) 
l. 


e2  ~ 


~_2_ 

d 


=  displacement  of  the  second  fringe  from  its 
initial  position  (inches) 

Once  the  temperatures  and  Tg ,  at  distances  y^  and 
y2  from  the  plate,  are  known,  then  a  second-order  curve 
can  be  passed  through  these  two  points: 

T  =  a  +  by  +  cy2 

where  a  =  T, 


T2  '  T1 

b  -  y  -  cyx 


and 


y2(T2  -  Tx)  -  y;L(T3  -  T.^ 
°  =  yly2^yl  "  y2^ 


The  local  heat  transfer  at  the  wall,  q  ,  is  then 


■  3T 

_  3 \r 


w  ay'y=o 


=  -  k  b 
w 


where  kw  =  thermal  conductivity  of  air  at  the  wall  tem¬ 
perature  (fiFArcg) 

The  total  heat  transfer,  Q,  is  determined  by  dividing  the 
plate  into  small  sections,  determining  the  local  heat 


transfer  for  each  section  and  numerically  integrating  the 
local  heat  transfers  over  the  plate. 

The  average  heat  transfer  is  determined  by 

h  -  Q 

■  a(th  -  v 

where  A  =  area  of  the  plate  (fta) 

T^  =  temperature  of  the  hot  plate  (°F) 

=  temperature  of  the  cold  plate  (°F) 

The  average  Nusselt  number,  Nu,  for  the  plate  is  then 


where  L  =  width  of  the  air  layer  between  the  two  plates 
(ft) 

k  =  thermal  conductivity  of  air  at  the  average 

4-  T  11 

H  C  ,  Btu  . 
temperature  - ^ -  ^hr~ft  °F^ 

The  above-mentioned  technique  for  determining  the  average 
Nusselt  number  has  been  programmed  for  data  processing, 
see  Appendix  B.  The  procedure  for  determining  local  heat 
transfer  at  the  cold  plate  is  similar. 

A  more  accurate  value  of  the  average  Nusselt  number 
can  be  determined  if  the  plate  is  subdivided  into  a  large 
number  of  sections  and  the  local  heat  transfer  determined 
for  each  of  those  sections.  For  the  four  geometries  having 
a  height  of  5.9"  or  more,  the  plate  was  divided  into  20 
sections,  and  for  the  geometry  having  a  2.786"  height,  the 
plate  was  divided  into  8  sections. 


The  most  difficult  task  in  the  whole  project  was  found 
to  be  the  measuring  of  the  distances  on  the  interferograms , 
which  were  unfortunately  not  very  clear.  There  were  40 
interferograms  to  be  analyzed  and  a  total  of  3712  measure¬ 
ments  to  be  made  on  those  interferograms.  Two  techniques 
were  tried  to  ease  the  problem  but  they  failed.  The  first 
technique  was  the  use  of  "densitometer"  which  is  available 
in  the  Aero-Propulsion  Laboratory  of  Wright-Patterson  U.S. 
Air  Force  Base.  The  densitometer  scans  a  negative  film  or 
a  transparency  for  density  variations  with  the  help  of  a 
thin  beam  of  light.  The  density  distribution  is  automat¬ 
ically  plotted  on  a  graph  with  a  magnified  distance  scale. 
The  method  did  not  work  because  if  too  thin  a  beam  of  light 
was  used  to  scan  the  negative,  then  the  densitometer  pro¬ 
vided  so  much  detail  of  the  density  variations  that  it 
became  impossible  to  find  the  center  of  a  fringe.  In  the 
case  of  a  thicker  beam  of  light,  the  densitometer  would 
not  detect  the  clustered,  thin  fringes  near  the  walls  of 
the  cavity.  The  second  technique  was  to  project  an 
enlarged  image  of  the  negative  on  a  screen  and  make  the 
measurements  with  a  scale.  This  technique  gave  very  bad 
results  and  was  also  abandoned.  The  measurements  were 
finally  made  with  the  help  of  a  microcomparator. 


IV.  RESULTS 


Interf erograms 

As  mentioned  earlier,  in  Section  III,  the  study  was 
carried  out  for  five  different  aspect  ratios.  Each  aspect 
ratio  was  studied  at  four  different  angles  of  inclination 
and  two  plate  temperature  differences.  The  results 
achieved  from  these  studies  are  included  in  this  section. 
Figures  10-15  include  the  interf erograms  for  the  different 
geometries  at  various  angles.  These  interferograms  are 
for  the  lower  temperature  difference  between  the  plates. 

The  rest  of  the  interferograms  are  in  the  custody  of  the 
Department  of  Aeronautics  and  Astronautics  at  AFIT. 

Figure  10  consists  of  the  interferograms  for  the 
5.93"  x  0.525"  geometry  (Aspect  Ratio  =  11.295).  Figure 
10(a)  shows  the  initial  finite  fringe  pattern  when  there 
is  no  temperature  difference  between  the  vertical  plates. 
The  fringes  have  to  be  adjusted  in  a  similar  way,  per¬ 
pendicular  to  the  plates,  whenever  a  cavity  is  placed  at 
a  new  angle.  Figure  10(b)  is  the  interf erogram  of  the 
vertical,  5-93"  x  0.525"  geometry  with  a  temperature  dif¬ 
ference  between  the  plates  of  24.6°F.  It  is  clear  from 
the  interferogram  that  the  fringes  in  the  central  part 
of  the  cavity  are  almost  straight  lines,  indicating  a 
linear  temperature  distribution  between  the  plates.  This 
means  that  in  most  of  the  central  portion  of  the  cavity, 
the  mode  of  the  heat  transport  is  predominantly  conduction. 


Some  convective  heat  transport  is  taking  place  at  both 
ends  of  the  cavity.  The  fringes  are  close- together  near 
the  bottom  of  the  hot  plate  and  the  top  of  the  cold  plate, 
indicating  the  presence  of  boundary  layers.  The  fluid 
flow  pattern  can  be  pictured  as  one  large  cell  with  the 
fluid  rising  near  the  hot  plate,  falling  near  the  cold 
plate.  Because  the  two  plate  boundary  layers  interfer, 
the  streamlines  are  nearly  parallel  to  the  bounding  sur¬ 
faces  everywhere  except  at  the  extreme  ends  of  the  cavity 
where  the  fluid  turns.  Thus,  the  heat  transfer  is  mostly 
conductive  except  at  the  ends. 

Figures  10(c),  (d),  and  (e)  are  the  interf erograms  of 
the  same  geometry  (5*93"  x  0.525")  at  60°,  30°,  and  0° 
angles  of  inclination.  All  of  these  interf erograms  can 
be  seen  only  partially.  To  find  the  average  heat  flux 
over  any  of  the  plates,  it  is  necessary  that  the  distri¬ 
bution  of  local  heat  fluxes  be  known  over  the  entire  plate. 
To  overcome  this  problem,  it  has  been  assumed  that  the 
local  heat  flux  distribution  over  the  top  half  of  the  hot 
plate  is  equivalent  to  that  on  the  lower  half  of  the  cold 
plate.  This  assumption  has  been  found  to  be  valid  by 
previous  studies  and  to  some  extent  by  the  present  study 
for  all  the  geometries  at  90° »  when  the  whole  cavity  is 
visible  on  the  interf erograms .  It  is  clear  from  Figure 
10(c)  that  the  mode  of  the  heat  transport  over  most  of 
the  central  part  of  the  cavity  is  still  conduction.  Fig¬ 
ures  10(d)  and  (e)  show  that  convection  is  the  predominant 
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means  of  heat  transport  over  the  entire  cavity,  even  though 
the  &T  for  the  0°-case  (Figure  10(e))  is  only  19.34?°F  as 
compared  to  24.5?5°F  of  the  90°-case  (Figure  10(b)). 

Figure  11  consists  of  the  interferograms  of  2.?86" 
x  3.00"  geometry  (Aspect  Ratio  =  0.929).  It  is  clear  from 
the  interferograms  that  in  all  of  the  cases,  it  is  con¬ 
vection  which  is  the  predominant  mode  of  heat  transfer. 

It  can  also  be  observed  that  in  all  the  cases,  except 
when  the  cavity  is  at  0°  (Figure  11(d)),  the  fringes  are 
very  close  together  near  the  bottom  of  the  hot  plate,  and 
they  separate  upwards  along  the  hot  plate.  The  same  phe¬ 
nomenon  can  be  observed  along  the  cold  plate  also,  but  in 
the  opposite  direction.  This  shows  that  a  boundary  layer 
is  formed  on  both  the  plates,  which  is  the  thinnest  at 
the  bottom  of  the  hot  plate  and  the  top  of  the  cold  plate. 
It  can  also  be  observed  from  all  the  interferograms  that 
the  fringes  in  the  core-region  between  the  two  boundary 
layers  are  almost  perpendicular  to  the  plates,  thus  having 
very  little  displacement  from  their  initial  position.  This 
means  that  there  is  very  little  conduction  in  the  core¬ 
region.  Since  the  heat  flux  is  concentrated  in  the  bound¬ 
ary  layer  regions  and  the  boundary  layers  are  the  thinnest 
at  the  bottom  of  the  hot  plate  and  the  top  of  the  cold 
plate,  it  can  be  predicted  from  the  interferograms  that 
the  rate  of  heat  flux  will  be  very  high  near  the  bottom 
of  the  hot  plate  and  the  top  of  the  cold  plate.  There  is 
some  strange  phenomenon  occurring  in  the  top  corner  of  the 
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core-region,  near  the  hot  plate,  and  in  the  bottom  corner 
of  the  core-region,  near  the  cold  plate,  which  is  as  yet 
unexplained.  It  can  be  observed  from  Figure  11(d)  that 
in  the  0°-case,  the  boundary  layers  seem  to  be  of  almost 
uniform  thickness  over  the  plates.  These  are  essentially 
stagnant  conduction  layers  separated  by  the  core-region 
where  energy  is  transported  by  convection. 

Figure  12(a)  consists  of  a  2.786"  x  3-00"  geometry 
placed  at  90°  but  with  a  diagonal  partition  inserted  into 
the  cavity.  The  partition  was  made  of  0.1"  thick  hard 
cardboard  with  both  sides  covered  with  white  paper.  The 
partition  was  placed  in  such  a  way  that  the  hot  plate 
remained  below  it.  By  the  insertion  of  the  partition, 
two  new  boundary  layers  have  been  introduced  into  the 
flow,  one  along  each  surface  of  the  partition.  The  tri¬ 
angular  cavities  formed  by  the  partition  have  more  stagnant 
areas  in  the  corners,  resulting  in  less  convection.  Both 
of  these  effects  decrease  the  heat  flux  between  the  plates. 

Figure  12(b)  is  the  interferogram  for  5-9"  x  3.0" 
geometry  placed  at  90°  and  with  a  aT  =  22.51°F.  It  is 
clearly  noticeable  that  the  fringe  pattern  seems  to  be 
similar  to  that  of  a  2.786"  x  3.00"  geometry;  therefore, 
the  rest  of  the  interf erograms  have  not  been  exhibited  in 
the  report. 

Figure  13  consists  of  the  interferograms  of  5* 9"  x 
1.5"  geometry  for  all  four  angles.  One  can  notice  from 
these  pictures  that  the  flow  is  definitely  of  the  boundary 


Figure  13 

a)  5-9"  x  1.5"  at  90°  c)  5-9"  x  1.5"  at  30° 

with  AT  =  2 5. 788° P  with  at  =  13.^35°F 

b)  5.9"  x  1.5"  at  60°  cl)  5-9"  x  1.5"  at  0° 

with  AT  =  20.522°F  with  AT  =  21.8l8°F 
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layer  type.  It  can  also  be  noticed  that  the  deflection 
of  the  fringes  in  the  core-region  becomes  less  as  one 
decreases  the  angle  of  inclination  from  90°  to  0°.  It 
can  probably  be  deduced  that  the  core  becomes  more  and 
more  isothermal  as  one  decreases  the  angle  of  inclination 
of  the  cavity. 

Figure  14(a)  is  the  interf erogram  of  5*9”  x  0.75" 
geometry  at  90°  and  with  aT  =  19.3^8^,  while  Figure  14(b) 
is  the  interferogram  of  the  same  geometry  at  the  same  angle 
but  with  at  =  33*8l8°F.  The  fringes  become  more  slanted 
and  closer  together  as  the  temperature  difference  between 
the  plates,  at,  is  increased.  Both  of  the  effects,  of 
course,  mean  higher  heat  flux.  Figures  14(c)  and  (d)  are 
the  interf erograms  at  two  different  aTs  for  the  60°-case 
of  the  same  geometry.  These  two  pictures  show  the  same 
effects  of  a  higher  at.  It  can  be  noticed  from  Figures 
14(a)  and  (c),  the  lower  aT  cases  for  both  the  angles, 
that  the  fringes  in  the  central  protion  of  the  cavity  cam 
be  approximated  as  straight  lines,  thus  showing  that  con¬ 
duction  may  be  a  significant  part  of  the  heat  transfer. 

While  there  seems  to  be  almost  no  core  region  for  the 
lower  at  cases  of  both  the  angles,  there  is  a  definite 
core-region,  which  is  almost  isothermal  for  the  higher 
aT  cases.  Figure  15  consists  of  interferograms  of  the 
same  geometry  at  30°  and  0°  with  lower  and  higher  aT. 

These  interferograms  show  that  the  flow  is  boundary  layer 
type  for  all  of  the  cases,  and  convection  is  the  predominant 


Figure  15 

a)  5.9"  X  0.75"  at  30°  c)  5-9"  x  0. 75"  at  0° 

with  AT  =  23.803°F  with  AT  =  l8.863°F 

b)  5. 91'  X  0.75"  at  30°  d)  5-9"  x  0.75"  at  0° 

with  AT  =  4l.608°F  with  AT  =  38.304°F 
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Local  Heat  Flux 

Local  heat  fluxes  were  calculated  at  twenty  positions 
along  the  hot  plate  and  the  cold  plate  for  the  geometries 
having  a  height  of  5*9"  and  at  eight  positions  for  the 
2.786"  x  3-00"  geometry.  These  positions  are  equally 
spaced  and  are  numbered  starting  from  the  bottom  of  the 
plate.  Local  heat  fluxes  for  two  vertical  geometries  are 
presented  in  Tables  I  and  II.  The  rest  of  the  local  heat 
flux  data  is  provided  in  Appendix  C. 


Table  I 

Local  Heat  Fluxes  for  5*93"  x  O.525"  Geometry 
Aspect  Ratio  =  11.295,  *T  =  24.5?5°F,  +  =  90° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
( Btu/hr  ft2) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2) 

1 

0.0 

10.3 

-  3-9 

2 

0.324 

26.9 

-  3.8 

3 

0.648 

20.6 

-  6.2 

4 

0.972 

16.5 

-  7-4 

5 

1.296 

13.8 

-  9.4 

6 

1.620 

12.9 

-  7.3 

7 

1.945 

9-7 

-  9.8 

8 

2.227 

10.1 

-10.2 

9 

2.593 

10.2 

-10.1 

10 

2.917 

10.1 

-10.3 

11 

3.014 

10.6 

-  9.8 

12 

3-338 

9.9 

-10.0 

13 

3.662 

9.2 

-  9.9 

14 

3.086 

9.6 

-10.0 

15 

4.310 

7.4 

-13.2 

16 

4.634 

9.0 

-13.9 

17 

4.958 

7.3 

-16.1 

18 

5.283 

6.6 

-21.3 

19 

5.607 

4.2 

-28.5 

20 

5.930 

4.3 

-10.1 
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Table  II 

Local  Heat  Fluxes  for  2.?86"  x  3.00"  Geometry 
Aspect  Ratio  =  O.929,  aT  =  24.l83°F,  t  =  90° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2) 

1 

0.0 

3.1 

-  3.8 

2 

0.398 

18.4 

-  4.2 

3 

0.796 

15-4 

-  8.8 

4 

1.194 

9-9 

-  9.4 

5 

1.592 

9.6 

-12.8 

6 

1.990 

6.4 

-17.5 

7 

2.388 

5.1 

-19.6 

8 

2.786 

5.3 

-  3.4 

Average  Heat  Flux  at  the  Hot  Plate  =  9.87 
Average  Heat  Flux  at  the  Cold  Plate  =  10. 72 


The  local  heat  flux  data  of  the  two  geometries  has 
been  presented  here  to  show  that  the  local  heat  fluxes  on 
correspondingly  opposite  positions  of  the  two  plates  are 
at  least  similar,  if  not  nearly  equal.  The  percentage 
difference  of  the  local  heat  fluxes  of  position  1  of  the 
hot  plate  and  position  20  of  the  cold  plate  of  the  5.93" 


x  0.525"  geometry  (see  Table  I),  for  example,  is  only  2.0 5$. 
The  main  concern  is  about  the  equality  of  the  local  heat 


fluxes  on  correspondingly  opposite  positions  at  the  top  of 
the  hot  plate  and  the  bottom  of  the  cold  plate,  because  in 
any  inclined  geometry  the  local  heat  fluxes  on  the  top  half 
of  the  hot  plate  have  been  substituted  for  the  local  heat 
fluxes  on  the  lower  half  of  the  cold  plate.  This  was,  of 
course,  not  necessary  for  the  2.?8 6"  x  3.00"  geometry, 
because  that  particular  geometry  was  completely  visible 
in  the  interferogram.  It  can  be  noticed  from  Table  I  that 
for  the  5*93"  x  0.525"  geometry  the  difference  between  the 
local  heat  fluxes  at  position  20  (top)  of  hot  plate  and 
position  1  (bottom)  of  the  cold  plate  is  11.5$.  The  dif¬ 
ference  for  the  same  two  corners  for  the  second  case  (see 
Table  II)  is  31-5^$  which  is  much  larger  than  the  other 
value.  Figure  1 6  shows  the  temperature  distributions  for 
a  2.786"  x  3. 00"  geometry,  placed  at  90°  and  with  a  T  = 
38.749°F.  It  is  clear  that  the  slope  of  the  temperature 
profile  at  the  top  of  the  hot  plate  is  more  than  the  slope 
at  the  bottom  of  the  cold  plate,  indicating  that  the  local 
heat  flux  at  the  top  of  the  hot  plate  is  higher  than  the 
heat  flux  at  the  bottom  of  the  cold  plate.  It  can  also  be 
noticed  from  Tables  I  and  II  that  there  is  a  lot  of  scatter 
in  the  data  of  local  heat  fluxes.  These  differences  are 
most  probably  due  to  random  errors  in  measuring  the  dis¬ 
tances  on  an  interferogram.  Significant  errors  are  also 
introduced  while  calculating  the  temperature  gradients, 
numerically.  In  addition,  complex  radiative  effects  may 
also  cause  some  difference  between  the  local  heat  fluxes 


at  the  top  of  the  hot  plate  and  the  bottom  of  the  cold 
plate.  Nevertheless,  the  local  heat  fluxes  at  the  top  of 
the  hot  plate  constitute  a  small  part  of  the  total  heat 
flux  from  the  hot  plate,  so  the  error  introduced  by 
replacing  these  be  the  local  heat  fluxes  at  the  bottom 
of  the  cold  plate  is  not  too  high. 

Figures  17-19  show  the  distribution  of  local  Nusselt 
numbers  on  the  hot  plate  for  5-9”  x  0.75”  geometry  (Aspect 
Ratio  =  7-867)  for  90°,  60°,  and  30°  angles  of  inclination, 
respectively. 


Local  Nusselt  Number  =  ^7*^ 

where  q  =  local  heat  flux  (Btu/hr  ft  ) 

L  =  thickness  of  air  layer  (ft) 

T  =  temperature  difference  between  the  plates  (°F) 
k  =  thermal  conductivity  of  air  (Btu/hr  ft  °F) 

It  can  be  observed  from  Figures  17-19  that  the  heat  flux 
at  the  very  bottom  of  the  hot  plate  is  low  and  it  rises 
immediately  to  a  peak  value  as  one  moves  a  small  distance 
up  the  hot  plate.  The  low  value  of  the  heat  flux  at  the 
very  bottom  of  the  hot  plate  is  due  to  a  stagnant  region 
in  the  corner  of  the  cavity.  The  heat  flux  rises  to  a 
high  value  as  one  goes  out  of  this  stagnant  region,  because 
at  the  bottom  of  the  hot  plate  the  boundary  layer  is  very 
thin.  The  local  heat  flux  then  keeps  decreasing  as  one 
moves  up  the  hot  plate  because  the  boundary  layer  develops 


•  .  ■ 


38 


STANCE  ALONG  THE  PLATE ( INCHES  ) 


ASPECT  RRTI0=  7.867 
PHI=  60. 


00  ’8  00'9  OO’V  00  *2 

aaauriN  iiassnN  iuoqi 


00*0 


Figure  18  .  Distribution  of  local  Nusselt  number  on  the 

hot  plate  for  5.9"  x  O.75* 
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and  becomes  thicker.  It  can  be  observed  from  these  plots 
that  as  the  angle  of  inclination  is  changed  from  90°  to  60° 
and  then  to  30°,  the  heat  flux  drops  off  more  steeply  after 
reaching  a  peak  value,  as  one  moves  up  the  hot  plate. 

After  the  drop,  the  heat  flux  becomes  more  uniform  along 
the  hot  plate.  The  reason  is  that  the  buoyant  force,  which 
always  works  in  the  vertical  direction  against  gravity,  has 
a  component  perpendicular  to  the  plate  when  the  angle  of 
inclination  is  any  other  than  90° •  This  component  of  the 
buoyant  force  pulls  the  air  in  the  boundary  layer  away  from 
the  plate,  thus  thickening  the  boundary  layer,  faster.  As 
the  angle  of  inclination  is  decreased  further,  this  effect 
increases  and  the  boundary  layer  develops  even  faster.  The 
faster  the  development  of  the  boundary  layer,  the  steeper 
would  be  the  drop  in  the  local  heat  flux.  As  the  boundary 
layers  on  the  hot  and  the  cold  plates  grow  faster,  they 
soon  join  each  other  and  the  heat  flux  becomes  almost  uni¬ 
form  over  both  the  plates. 

Figure  20  is  the  local  Nusselt  number  distribution  on 
the  hot  plate  for  5 •9"  x  0 .75"  geometry  (Aspect  Ratio  = 
7.867)  at  0°.  It  can  be  observed  from  this  plot  that  the 
heat  flux  is  low  at  the  very  bottom  of  the  hot  plate  due  to 
the  stagnant  region,  it  rises  to  a  high  value  as  one  goes 
up  the  hot  plate,  remains  almost  uniform  over  most  of  the 
plate  and  then  drops  down  as  one  approaches  the  stagnant 
region  at  the  top  of  the  plate.  The  exact  nature  of  the 
flow  in  a  horizontal  cavity  is  quite  complex.  It  is  quite 


evident  from  the  interferograms  and  Figure  20  that  there  is 
almost  uniform  heat  flux  over  the  plates  and  most  of  the 
temperature  change  takes  place  in  a  thin  layer  near  the 
plates,  while  the  core  is  almost  isothermal.  Figure  21  is 
the  local  Nusselt  number  distribution  over  the  cold  plate 
for  the  same  geometry  at  90°.  As  can  be  seen,  it  is  exactly 
opposite  to  Figure  1?. 
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Average  Heat  Flux 


The  average  heat  fluxes  are  provided  in  Tables  III  - 
VII  in  the  form  of  average  Nusselt  numbers,  Nu,  where  the 
average  Nusselt  number  Nu  is  based  on  L,  the  width  of  the 
air  layer.  The  average  Nusselt  numbers  have  been  presented 
against  either  Grashof  numbers  in  the  case  of  90°  inclina¬ 
tion  (Tables  III  and  IV)  or  Rayleigh  numbers  in  the  case  of 
60°,  30°  and  0°  inclinations  (Tables  V,  VI,  and  VII).  Both 
the  Grashof  numbers  and  the  Rayleigh  numbers  are  based  on 
L,  the  width  of  the  air  layer.  The  theoretical  Nusselt 
numbers  given  in  Table  III,  which  consists  of  low  Grashof 
numbers  (conductive  regime),  were  calculated  by  the  empir¬ 
ical  equation  (1). 

Nu  =  1  +  0.00166  ^  (Gr)0,9  (1) 


Table  IV,  which  consists  of  high  Grashof  numbers  (boundary 
layer  regime),  provides  theoretical  Nusselt  numbers  which 
were  calculated  by  the  empirical  equation  (2). 

Nu  =  0.119  (Gr)0,3  (^)0,1  (2) 


Equations  (1)  and  (2)  were  extracted  from  Reference  8.  The 
theoretical  Nusselt  numbers  provided  in  Tables  V  -  VII  were 
calculated  by  the  empirical  equation  (3). 


Nu  =  1  +  1.44  [1  - 


1708 


Ra  cos  $ 


1*  (sin  1.8$) 

3  Do  oo 


1.6 


Ra  cos  <f 


1708, 


+  f Ra  cos  <Kl/3 

' 555o '  11 

where  [ ]  is  defined  by: 


(3) 
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tJtAS 


[x3  =  ( Ixl  +  x)/2 

Equation  (3)  was  extracted  from  Reference  6 


Table  III 

Average  Nusselt  Numbers  for 
Low  Grashof  Numbers,  $  =  90° 


Grashof 

Number 

Nu 

(Experimental) 

Nu 

(Theoretical) 

% 

Difference 

4613.0? 

1.34 

1.29 

3.8 

6229 • 54 

1.27 

1.38 

8.3 

17497.02 

2.21 

2.39 

7.8 

Table  IV 

Average  Nusselt  Numbers  for 

High  Grashof  Numbers,  $  =  90° 

Grashof 

Nu 

Nu 

% 

Number 

(Experimental) 

(Theoretical) 

Difference 

107662.20 

3.34 

3-36 

0.6 

152281.42 

3.90 

3.73 

5.2 

766760.97 

6 . 63 

6.48 

2.3 

843023.22 

7.14 

7.19 

0.7 

1244511.49 

7.83 

7.49 

4.4 

1287084.44 

8.68 

8.16 

6.2 
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Table  V 


Average  Nusselt  Numbers  for  $  =  60° 


Rayleigh 

Number 

Nu 

(Experimental ) 

Nu 

(Theoretical) 

% 

Difference 

2994.1 

0.99 

1.00 

1.0 

4386.4 

1.22 

1.09 

11.2 

8403.7 

1.51 

1.53 

1.3 

13563.2 

1.70 

1.88 

10.1 

62647.2 

2.74 

3.04 

10.4 

111669.0 

3.55 

3.48 

2.0 

455359.8 

5.27 

4.81 

9.1 

630896.8 

5.57 

5.21 

6.7 

1099864.0 

5.74 

5.98 

4.2 

Table  VI 

Average  Nusselt  Numbers  for  $  =  30° 

Rayleigh 

Nu 

Nu 

% 

Number 

(Experimental) 

(Theoretical) 

Difference 

3496.3 

1.23 

1.37 

10.8 

4963.1 

1.50 

1.62 

7.7 

9016.9 

1.90 

2.05 

7.6 

43454.9 

3.12 

3.19 

2.2 

113075.0 

3.46 

3.96 

13-5 

605935.5 

6.44 

5.91 

8.6 

877750.9 

6.08 

6.51 

6.8 

1150561.1 

6.89 

6.99 

1.4 
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Table  VII 


Average  Nusselt  Numbers  for  *  =  0° 


Rayleigh 

Number 

Nu 

(Experimental ) 

Nu 

(Theoretical) 

% 

Difference 

2587.3 

1.47 

1.49 

1.4 

4746.8 

1-95 

1.92 

1.6 

7201.0 

1-95 

2.17 

10.7 

67578.2 

3-30 

3.67 

10.6 

560077.7 

6.44 

6.02 

6.7 

609960.7 

6 . 60 

6.15 

7-1 

829539-9 

6.16 

6.66 

7.8 

936604.7 

6.74 

6.87 

1.9 

The  results  presented  in  Tables  III  -  VII  are  given  in 
Figures  22-26,  respectively.  Experimental  and  theoretical 
Nusselt  numbers  are  plotted  against  Grashof  numbers  in  the 
case  of  90°  (Figures  22  and  23),  and  these  have  been  plotted 
against  Rayleigh  number  x  cos  $  for  60°,  30°,  and  0°  angles 
of  inclinations.  It  can  be  seen  that  the  three  empirical 
equations  are  within  11%  difference  from  the  experimental 
results .  Considering  the  random  errors  in  the  experimental 
results,  this  error  can  be  accepted  and  the  empirical  equa¬ 
tions  can  be  expected  to  give  good  results  for  Grashof 
numbers  up  to  1  million. 

Figure  27  shows  the  dependence  of  average  Nusselt  num¬ 
ber  on  the  angle  of  inclination,  ♦ ,  for  a  Grashof  number  of 
200,000  and  Aspect  Ratio  of  3-933  (5-9”  x  1.5"  geometry). 

The  Nusselt  number  does  not  vary  much  as  the  angle  of  incli¬ 
nation,  ♦,  is  reduced  from  90°  to  60°.  However,  there  is  a 


- 


definite  rise  of  about  12%  in  the  Nusselt  number  as  $  is 
reduced  from  60°  to  30° >  and  there  is  again  very  little 
change  as  $  is  further  reduced  from  30°  to  0°.  Probably 
the  flow  mechanism  in  the  cavity  does  not  alter  much  if 
the  angle  of  inclination  of  the  cavity  is  changed  from 
90°  to  60°.  There  seems  to  be  a  transition  from  the 
boundary  layer  type  flow  of  the  vertical  geometry  to  the 
local  cells  type  flow  of  the  horizontal  geometry,  as  the 
angle  of  inclination  is  reduced  from  60°  to  30°.  The 
changes  in  the  flow  seem  to  be  not  much  again,  if  the 
angle  of  inclination  is  further  reduced  from  30°  to  0°. 

It  can  be,  however,  noticed  that  changes  in  the  angle  of 
inclination  of  an  average  aspect  ratio  geometry  do  not 
bring  about  drastic  changes  in  the  average  Nusselt  number, 
if  the  Grashof  number  is  kept  constant. 
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heoretical  and  experimental  Nusselt  numbers 


Figure  25  .  A  plot  of  theoretical  and  experimental  Kusselt  numbers 

vs  (Rayleigh  number  x  cos  4> )  . 
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For  the  partitioned  case 


Grashof  number  =  1267224.8 

a 

Average  Nusselt  number  =  3.65 

while  for  the  case  without  the  partition  (same  geometry  and 
<t> ) 


Grashof  number  =  1287084.4 

n  rrn  T\7  7  t  o  1  +■  _  ft  Aft 


V.  Conclusions 


Based  on  the  results  of  this  study,  the  following  cun 
elusions  are  drawn: 

1.  The  average  heat  flux  across  a  rectangular  cavity 
whether  vertical  or  inclined,  is  primarily  a  function  of 
the  Grashof  number  or  Rayleigh  number. 

2.  For  a  constant  Grashof  number,  the  average  heat 
flux  across  a  rectangular  cavity  does  not  change  apprecia¬ 
bly  by  either  changing  the  aspect  ratio  or  the  angle  of 
inclination. 

3.  The  correlation  equations  developed  by  Eckert  for 
the  evaluation  of  average  heat  flux  across  a  vertical,  rec 
tangular  cavity  and  the  correlation  equation  developed  by 
Hollands  to  determine  the  average  heat  flux  across  an 
inclined,  rectangular  cavity,  approximate  the  experimental 
data  within  +  11$,  up  to  a  Grashof  number  of  1.2  million. 

4.  The  introduction  of  a  diagonal  partition  into  a 
rectangular  cavity  decreases  the  heat  flux  more  than  half. 


VI.  Recommendations 


It  is  recommended  that: 

1.  Further  experimental  study  should  be  carried  out 
on  a  partitioned  cavity,  because  the  introduction  of  a 
partition  decreases  the  heat  flux  significantly. 


2.  A  future  study  of  heat  flux  should  be  carried  out 
with  the  help  of  a  Moire  schlieren  apparatus,  which  might 
provide  a  clearer  fringe  pattern. 
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APPENDIX  A 


Finite  Fringe  Method 

Two  beams  of  light,  having  sinusoidal  wave  patterns 
with  the  same  wavelength  ,  when  combined,  cancel  each  other 
if  the  two  are  out  of  phase  by  l80°  and  complement  each 
other  if  in  phase.  Therefore,  a  series  of  coherent  beams 
parallel  to  and  in  phase  with  each  other,  joining  with 
another  coherent  series  of  beams,  will  produce  either  an 
area  of  light,  or  an  area  of  dark,  on  a  surface  in  their 
path,  depending  on  whether  the  wave  patterns  of  the  two 
beams  are  in  phase  or  out  of  phase  by  l8o°  when  they  join. 

In  the  finite  fringe  mode  of  operation,  mirrors, 
and  Mg ,  and  splitter  plates,  S.^  and  Sg,  (refer  to  Figure  7) 
are  set  at  approximately  45°  angles.  Then,  minor  adjust¬ 
ments  are  made  to  first  produce  an  infinite  fringe  pattern 
having  equal  path  lengths  for  the  test  section  and  refer¬ 
ence  beam.  The  splitter  plate,  Sg ,  is  then  rotated  through 
a  small  angle,  e,  as  shown  in  Figure  28. 
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Test  Beam 


Figure  28.  Test  and  Reference  Beams 
in  Finite  Fringe  Mode. 

The  first  interference  fringe,  from  the  axis  of  rota¬ 
tion,  forms  where  the  test  section  beam  path  length  has 
been  shortened  by  x0/2,  and  the  reference  beam  path  length 
has  been  lengthened  by  xo/2,  where  is  the  wavelength  of 
the  reference  beam  of  light.  From  the  geometry  in  Figure 
28 

cos  45° 

e  -  - ; - - - ; - 

d  x  1.  _  o 

- -  +  _ o  cos  45 

cos  45  2 

Now  d  >>  A  therefore  the  second  term  in  the  denominator 


can  be  neglected. 


e 


( A-l ) 


A 

O  I,  r-  o  , 

—  cos  ^5  *Q 

d  _ 

cos  45° 

An  optical  amplification  factor,  F,  is  defined  as  the 
ratio  of  the  distance  between  fringes,  d,  in  a  finite 
fringe  pattern,  and  the  difference  in  beam  path  lengths 
between  fringes,  Aq. 


o 


Since  Aq  =  Cqt,  the  distance  between  finite  fringes,  d,  can 
be  expressed  in  terms  of  the  time  difference  between  the 
two  beams,  which  is  the  period  of  the  light  beam,  -r . 

d  =  FCqt  ( A-3) 

Under  test  conditions,  the  difference  in  time  for  the  test 
beam  and  the  reference  beam  is  changed  by  the  change  in  the 
velocity  of  light  in  the  test  section. 

AT  =  CT  "  C  (A-4) 

where  L  is  the  length  of  the  test  section.  From  equation 
(A-3),  this  time  difference  produces  a  change  in  fringe 
position,  Ad 

Ad  =  FCqAt  ( A-5) 

Dividing  Ad,  equation  (A-5),  by  d,  equation  (A-3);  and 
expressing  the  time  change  as  a  multiple,  e,  of  the  period 
of  radiation 


( A-6) 
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W  f 


t 


From  equation  (A-4)  and  equation  (A-6) 

T  T  A 

L  L  vac 

C  _C~eT_eC 
o  vac 


(A-7) 


Rearranging  and  introducing  the  definition  of  index  of 
C. 

refraction,  n  = 


vac 


(n0  -  n)L  -  elvac 


(A-8) 

The  index  of  refraction  is  related  to  gas  density  by  the 
empirical  Dale-Gladstone  constant,  K^. 


(nQ  -  n)  =  -pKdq 

Combining  equation  (A-8)  and  equation  (A-9) 

c  x 

vac 

P  =  '  KDGL 


(A-9) 


( A-10) 


The  difference  in  gas  density  between  two  positions  in  the 
test  section,  corresponding  to  two  points  on  an  interfer¬ 
ence  fringe,  depends  on  ie  between  them.  Let 


p  1  —  —  e 


vac 


vac 


1  K^L  and  p2  (el  +  Ae)  K_.r,L 


( A-ll ) 


DG 


DG 


Then 


vac 

KDGL 


For  an  ideal  gas  and  constant  pressure 

^  X 


( A-12 ) 


RT, 


RT„ 


=  -A  e 


vac 

KDGL 


"2  “*1 

Rearranging  equation  (A-13)  one  gets 

T 


(A-13) 


vac,  m 

PK^L  ti 


(A-l4) 


^*■111  ■!  I  »•  m  fr.  Ai?!  *  r. 


APPENDIX  B 


Computer  Program  for  Data  Reduction 
This  program  is  for  the  calculation  of  the  total  heat 
transfer  and  local  heat  transfers  at  20  sections  of  the  hot 
plate  for  any  of  the  4  geometries  having  a  height  of 
5-9"  or  more.  For  the  calculation  of  local  and  total  heat 
transfers  at  the  cold  plate,  the  only  change  required  is 
to  change  the  formulae  for  the  calculation  of  temperatures 
at  the  first  and  second  fringes.  In  that  case: 


T2  =  i  -  «lClTl 


and 

T1 

T3  ~  i  -  e2c1t1 

where  T1  =  temperature  of  the  cold  plate  (°R) 

For  the  geometry  having  2.786"  height,  only  eight 
sections  were  used  to  numerically  integrate  the  local  heat 
transfers  over  the  plates  to  calculate  the  total  heat 
transfer . 


Input  The  program  requires  the  following  parameters  as 
input : 

D  =  width  of  air  layer  (inches) 

T^  =  temperature  of  the  hot  plate  (°R) 

P  =  atmospheric  pressure  (pounds/ft2) 


K  =  thermal  conductivity  of  air  at  the  wall  tem¬ 
perature  (Btu/hr  ft  °F) 


CENT  =  width  of  the  central  section  on  the  plate 
Note :  The  plate  has  been  divided  into  20  sections 


for  the  purpose  of  numerical  integration  of  local  heat 
transfers  on  these  sections  to  get  total  heat  transfers. 
The  height  of  these  sections  is  as  follows: 

Sections  1  and  20  =  0.1  inches  on  picture 
Sections  2-9  and  12-19  =  0.2  inches  on  picture 
Sections  10  and  11  =  "cent"  inches  on  picture 


where 

CENT 

FF 


height  of  plate  on  picture-[2(0.1)+l6(0.2) ]  inches 
magnification  factor 

_ actual  height  of  the  air  layer _ _ 

height  of  the  air  layer  on  interf erogram 


L15L2  =  actual  displacements  of  the  first  and  the  second 
fringes  (inches) 

=  displacements  on  the  interf erogram  x  FF 
=  actual  distances  of  the  first  and  the  second 
fringes  from  the  wall  (feet) 


A  listing  of  the  program  is  as  follows: 


PROGRAM  TH3  (INPUT, OUTPUT) 

REAL  K 
REAL  LI 
REAL  L2 

20  READ*, D,T1,P,K, CENT, FF 
IF(D.EQ.0)G0  TO  25 
N=1 

5  READ* , LI , L2 , Y1 , Y2 

PRINT* , "Ll=" , LI , "L2=" ,L2 , "Yl=" , Y1 , "Y2=" , Y2 
EPS1=L1/D 
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C1=0. 026195/p 
T2=Tl/( 1 . 0+EPSl*Cl*Tl ) 

EPS2=L2/D 

T3=T1/ ( 1 . 0+EPS2*Cl*Tl ) 

C= ( Y2*( T2-T1 ) -Yl*( T3-T1 ) ) /( ( Y1-Y2 ) *Y1*Y2 ) 

B= ( ( T2-T1 ) /Y1 ) -( C*Y1 ) 

QLOC=-K*B 

PRINT* , "LOCAL  HEAT  TRANSFER  AT  HOT  PLATE=",QLOC 

IF( N . GT . 1 .AND . N. LT . 10 ) GO  TO  10 

IF( N . GE. 10 . AND . N . LT . 12 ) GO  TO  15 

IF( N. GE . 12 . AND . N . LT . 20) GO  TO  10 

IF( N . GE. 20 ) GO  TO  30 

QT0T=FF*0 . 1*QL0C/12 . 0 

N=N+1 

GO  TO  5 

10  QT0T=QT0T+0.2*FF*QL0C/12.0 

N=N+1 
GO  TO  5 

15  QT0T=QT0T+FF*CENT*QL0C/12 . 0 
N=N+1 
GO  TO  5 

30  QT0T=QT0T+FF*0 . 1*QL0C/12 . 0 

PRINT*, "TOTAL  HEAT  TRANSFER  AT  HOT  PLATE=" , QTOT 
GO  TO  20 

25  STOP"END  OF  PROGRAM" 

END 
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Appendix  C 


Local  Heat  Fluxes 


1 .  Angle  of  Tncli  nation,  $  =  90° 


Table  VIII 

Local  Heat  Fluxes  for  b-93"  x  0.525"  Geometry 
Aspect  Ratio  =  11.295,  AT  =  34.5l4°F,  ---  905 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 
Hot  Plate 
( Btu/hr  ft2 ) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2 ) 

1 

0.0 

11.8 

-  6.9 

2 

0.324 

25.6 

-  4.6 

3 

0.648 

23.8 

-  6.5 

4 

0.972 

22.0 

-  7.2 

5 

1.296 

21.2 

-15.6 

6 

1.62 

17.3 

-15.3 

7 

1.945 

10.7 

-17.0 

8 

2.269 

15.7 

-12.9 

9 

2.593 

16.1 

-17.0 

10 

2.917 

10.6 

-  9.2 

11 

3.014 

10.6 

-  9-2 

12 

3.338 

19.2 

-14.2 

13 

3.662 

l4 .6 

-]4.2 

l4 

3.986 

16.2 

-32.3 

15 

4.310 

9.8 

-15.2 

16 

4.634 

13.3 

-20.7 

17 

4.958 

7.4 

-21.4 

18 

5.283 

6.4 

-23.1 

19 

5.607 

5.1 

-23.4 

20 

5.930 

7-2 

-11.3 

Total  Heat 

Transfer  Rate 

at  Hot  Plate  = 

7.09  Btu/hr 

Total  Heat 

Transfer  Rate 

at  Cold  Plate 

=  -7.48  Btu/hr 
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Table  IX 


Local  Heat  Fluxes  for  2. 786”  x  3-00"  Geometry 
Aspect  Ratio  =  O.929,  at  =  38.844°F,  *  =  90° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2) 

1 

0.0 

11.3 

-  7.7 

2 

0.398 

33.9 

-14.1 

3 

0.796 

33.4 

-15.9 

4 

1.194 

17.7 

-18.2 

5 

1.592 

20.8 

-15.7 

6 

1.990 

16 . 6 

-28.6 

7 

2.388 

16.2 

-32.6 

8 

2.786 

10.8 

-  9.2 

Total.  Heat  Transfer  Rate  at  Hot  Plate  -  4.96  Btu/hr 
Total  Heat  Transfer  Rate  at  Cold  Plate  -  -4.43  Btu/hr 


Table  X 


Local  Heat  Fluxes  for  5-9"  x  3.00"  Geometry 
Aspect  Ratio=1.96?,  aT1=22.51°F,  aT2=38.94°F,  4>=90° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aT, 
(Btu/hr  ft2  ) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aTj 
(Btu/hr  ft2  ) 

1 

0.0 

4.5 

22.5 

2 

0.314 

9.4 

33.5 

3 

0.628 

12.4 

18. 7 

4 

0.942 

12.7 

19.1 

5 

1.257 

12.5 

33-9 

6 

1-571 

10.9 

25.6 

7 

1.885 

10.3 

18.6 

8 

2.199 

10.0 

22.1 

9 

2.513 

11.4 

16.0 

10 

2.827 

7.2 

23.9 

11 

3.072 

12.7 

22.0 

12 

3-386 

10.9 

I8.9 

13 

3.701 

9-2 

20.8 

14 

4.015 

7.0 

16.7 

15 

4.329 

8  „9 

10.5 

16 

4. 643 

8.3 

14.4 

17 

4.957 

5-0 

10.7 

18 

5.271 

5.1 

10.4 

19 

5.585 

2.3 

3.8 

20 

5.900 

3-7 

2.4 

Total  Heat  Transfer 
Total  Heat  Transfer 


Rate  at  Hot  Plate  for  AT-^ 

Rate  at  Hot  Plate  for  AT^ 

2 


4.43  Btu/hr 
9.18  Btu/hr 
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Table  XI 


Local  Heat  Fluxes  for  5*9"  x  1.5"  Geometry 
Aspect  Ratio  =  3-933,  aT  =  25.788°F,  <t»  =  90° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  ft2 ) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2 ) 

1 

0.0 

6.2 

-  2.5 

2 

0.318 

33.9 

-  1.0 

3 

0.635 

23.9 

-  1.0 

4 

0.952 

19.8 

-  3.2 

5 

1.270 

14.5 

-  3.4 

6 

1.587 

13.6 

-  4.1 

7 

1.905 

13.3 

-  4.1 

8 

2.222 

12.2 

-  6.1 

9 

2.540 

12.0 

-  8.4 

10 

2.857 

11.9 

-10.4 

11 

3.044 

11.8 

-11  .2 

12 

3-362 

8.6 

-10.7 

13 

3.679 

6.7 

-12.3 

14 

3-997 

4.3 

-  6.8 

15 

4.314 

3.5 

-12.2 

16 

4.632 

2.9 

-14.0 

17 

4.949 

2.4 

-19.9 

18 

5.267 

1.0 

-22.1 

19 

5.584 

1.2 

-32.7 

20 

5.900 

3.7 

-  5.9 

Total  Heat 

Transfer  Rate 

at  Hot  Plate  =  5 

.09  Btu/hr 

Total  Pleat 

Transfer  Rate 

at  Cold  Plate  =  ■ 

-4.89  Btu/hr 

Table  XII 


Local  Heat  Fluxes  for  5-9"  x  1.50"  Geometry 
Aspect  Ratio  =  3-933,  & T  =  37-505°F,  $  =  90° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 
Hot  Plate 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2) 

1 

0.0 

3.6 

-  4.2 

2 

0.318 

42.4 

-  2.9 

3 

0.635 

31.9 

-  3.5 

4 

0.952 

35-4 

-  6.7 

5 

1.270 

30.4 

-  8.2 

6 

1-587 

31.1 

-10.5 

7 

1.905 

25.6 

-  9.4 

8 

2.222 

19.1 

-11.4 

9 

2.540 

18.7 

-14.1 

10 

2.857 

17.5 

-14.7 

11 

3 . 044 

15.2 

-16.2 

12 

3.362 

15.1 

-20.0 

13 

3.680 

11.9 

-19.0 

14 

3-997 

8.9 

-24.1 

15 

4.314 

8.8 

-33.1 

16 

4.632 

8.1 

-26.9 

17 

4.949 

8.1 

-32.7 

18 

5.267 

3-7 

-30.0 

19 

5 . 584 

3.4 

-39.4 

20 

5.900 

5-2 

-  3.2 

Total  Heat 

Transfer  Rate 

at  Hot  Plate  = 

8.87  Btu/hr 

Total  Heat 

Transfer  Rate 

at  Cold  Plate  - 

-8.52  Btu/hr 
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Table  XIII 


Local  Heat  Fluxes  for  5*9"  x  0.75"  Geometry 
Aspect  Ratio  =  7.867,  AT  =  33.818°F,  <t>  =  90° 


Posi tion 
Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  ft2 ) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2 ) 

1 

0.0 

21.3 

-  3.4 

2 

0.315 

45.1 

-  2.4 

3 

0.630 

41.9 

-  2.6 

4 

0.944 

42.2 

-  3.4 

5 

1.259 

40.3 

-  5-1 

6 

1.574 

26.2 

-  5.3 

7 

1.889 

25.7 

-  5.9 

8 

2.204 

19.8 

-  7-8 

9 

2.518 

15.2 

-  9-8 

10 

2.833 

16.2 

-12.2 

11 

3.066 

12.8 

-14.4 

12 

3.381 

10.1 

-15.9 

13 

3.696 

8.0 

-20.0 

14 

4.011 

6.1 

-20.5 

15 

4.325 

5-6 

-2  6.7 

16 

4.640 

5.7 

-43.4 

17 

4.955 

2.3 

-47.8 

18 

5.270 

2.7 

-40.3 

19 

5-585 

3.1 

-52.1 

20 

5.900 

3.7 

-20.8 

Total  Heat  Transfer  Rate  at  Hot  Plate  =  8.86  Btu/hr 
Total  Heat  Transfer  Rate  at  Cold  Plate  =  -9.04  Btu/hr 


2 .  Angle  of  Inclination,  4>  -  60° 


Table  XIV 

Local  Heat  Fluxes  for  5 •93"  x  0.525"  Geometry 
Aspect  Ratio=11.295,  AT] =22 . 213°F ,  aT?=33 • 3?3°F ,  $=60° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  ATj 
(Btu/hr  ft5 ) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  at2 
(Btu/hr  ft2  ) 

1 

0.0 

7-6 

7.1 

2 

0.324 

11.5 

23.3 

3 

0.648 

12.4 

44.5 

4 

0.972 

10.1 

16.1 

5 

1.296 

10.0 

12.9 

6 

1.620 

8.9 

13.9 

7 

1.945 

6.9 

12.7 

8 

2.269 

5.4 

10.9 

9 

2.593 

5.9 

10.5 

10 

2.917 

4.5 

10.9 

11 

3.014 

5.7 

10.2 

12 

3-338 

6.5 

12.4 

13 

3.662 

7.4 

11.5 

14 

3-986 

7.0 

13.0 

15 

4.310 

7.6 

10.2 

16 

4.634 

8.1 

11.0 

17 

4.958 

7.0 

8.6 

18 

5.283 

6.8 

6 . 6 

19 

5.607 

3.4 

5-8 

20 

5.930 

3.7 

34.2 

Total  Heat  Transfer  Rate  at  Hot  Plate  for  AT-^ 
Total  Heat  Transfer  Rate  at  Hot  Plate  for  ATg 


3.64  Btu/hr 
6.84  Btu/hr 


75 


Table  XV 


Local  Heat  Fluxes  for  2.786"  x  3*00"  Geometry 
Aspect  Ratio  =  O.929,  aT  =  26.375°P,  ♦  =  60° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2) 

1 

0.0 

1.8 

-  0.7 

2 

0.398 

13.5 

-  3.4 

3 

0.796 

14.0 

-  8.1 

4 

1.194 

13.7 

-10.3 

5 

1.592 

9-3 

-14.0 

6 

1»990 

5.8 

“13.7 

7 

2.388 

4.3 

-13.7 

8 

2.786 

1.1 

-  1.7 

Total  Heat  '  ransfer  Rate  at  Hot  Plate  =  2.06  Btu/hr 
Total  Heat  Transfer  Rate  at  Cold  Plate  =  -2.34  Btu/hr 


Table  XVI 


Local  Heat  Fluxes  for  2.?8 6"  x  3*00"  Geometry 
Aspect  Ratio  =  O.929,  aT  =  49.859°F,  41  =  6 0° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2) 

1 

0.0 

0.9 

-  6.5 

2 

0.398 

23.9 

-12.0 

3 

0.796 

27.0 

-14.8 

4 

1.194 

20.8 

-18.5 

5 

1.592 

20.5 

-21.3 

6 

1.990 

13.8 

-25.6 

7 

2.388 

12.4 

-24.9 

8 

2.786 

7.0 

-  1.4 

Total  Heat  Transfer  Rate  at  Hot  Plate  -  4.06  Btu/hr 
Total  Heat  Transfer  Rate  at  Cold  Plate  =  -4.02  Btu/hr 


Table  XVII 


Local  Heat  Fluxes  for  5*9"  x  3-00  Geometry 
Aspect  Ratio  =  I.967,  aT  =  18.849°F,  4.  =  60° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  ft2  ) 

1 

0.0 

3-6 

2 

0.314 

3.7 

3 

0.628 

8.7 

4 

0 . 942 

9.3 

5 

1.257 

8.8 

6 

1.571 

7.3 

? 

1.885 

7.4 

8 

2.199 

6.4 

9 

2.513 

2.3 

10 

2.827 

4.4 

11 

3.072 

5-0 

12 

3.386 

7-3 

13 

3.701 

7.1 

14 

4.015 

7.9 

15 

4.329 

6.2 

16 

4.643 

7.3 

17 

4.957 

4.8 

18 

5.271 

3.8 

19 

5.585 

2.7 

20 

5.900 

1.6 

Total  He'at  Transfer  Rate  at  Hot  Plate  =  2.96  Btu/hr 


Table  XVIII 


Local  Heat  Fluxes  for  5 •9"  x  1 . 5”  Geometry 

Aspect  Rati0=3.933,  a'1’1=20 . 522°F ,  at2=38.566°F,  *=60° 


Position 

Number 

Distance  Along 
the  plate  "x” 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  AT, 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aTj 
(Btu/hr  ft2) 

1 

0.0 

2.0 

5-7 

2 

0.318 

15.2 

58.7 

3 

0.635 

13.8 

58.6 

4 

0.952 

13.8 

36.9 

5 

1.270 

6.8 

24.1 

6 

1.587 

7.5 

10.3 

7 

1.905 

5-3 

11.5 

8 

2.222 

8.3 

5-2 

9 

2.540 

8.4 

8.4 

10 

2.857 

6.8 

14.5 

11 

3.044 

5-2 

14.5 

12 

3.362 

5-5 

10.0 

13 

3.679 

5-5 

8.7 

14 

3.997 

6.1 

14.5 

15 

4.314 

3.6 

9-7 

16 

4.632 

2.6 

5.7 

1? 

4.949 

3.4 

5-8 

18 

5.267 

2.0 

5.9 

19 

5-584 

1.9 

1.9 

20 

5.900 

8.2 

8.1 

Table  XTX 


Local  Heat  Fluxes  for  5*9”  x  0.75"  Geometry 

Aspect  Ratio=?.867,  aTx=21 . 962°F ,  aT2=37.435°F,  $=60° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aT, 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aT2 
(Btu/hr  ft2) 

1 

0.0 

3.4 

0.0 

2 

0.315 

22.9 

71.3 

3 

0.630 

18.0 

35.8 

4 

0.944 

17.3 

27.9 

5 

1.259 

12.4 

21.2 

6 

1.574 

12.2 

15.5 

7 

I.889 

8.1 

15.2 

8 

2.204 

6.8 

14.0 

9 

2.518 

5-4 

15.1 

10 

2.833 

6.7 

14.4 

11 

3.066 

6.4 

11.7 

12 

3.381 

6.5 

11.5 

13 

3.696 

5.6 

8.2 

14 

4.011 

5.0 

8.4 

15 

4.325 

5.2 

8.2 

16 

4 . 640 

3.1 

3-2 

17 

4.955 

2.3 

3.0 

18 

5.270 

1.6 

3.4 

19 

5-585 

1.0 

1.2 

20 

'  5.900 

5.9 

9.0 

Total  Heat  Transfer  Rate  at  Hot  Plate  for  AT^  =  3-92  Btu/hr 
Total  Heat  Transfer  Rate  at  Hot  Plate  for  aT2  =  7.62  Btu/hr 
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3 •  Angle  of  Inclination,  »  =  30° 


Table  XX 

Local  Heat  Fluxes  for  5*93"  x  0.525"  Geometry 
Aspect  Ratio=11.295.  ^=26. 151°F,  aT2=39.00°F,  $=30° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aT, 
(Btu/hr  ft2 ) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aT2 
(Btu/hr  ft2) 

1 

0.0 

14.1 

25.9 

2 

0.324 

15.7 

32.3 

3 

0.648 

14.0 

20.9 

4 

0.972 

12.0 

17.5 

5 

1.296 

11.? 

16.1 

6 

1.620 

8.5 

17.3 

7 

1.945 

8.4 

16.2 

8 

2.269 

15.8 

20.8 

9 

2.593 

14.2 

25.9 

10 

2.917 

13.0 

28.4 

11 

3.014 

9.5 

20.8 

12 

3.338 

6.9 

16 . 8 

13 

3.662 

7.6 

18.4 

14 

3-986 

9-5 

17.7 

15 

4.310 

14.1 

18.3 

16 

4.-634 

11.3 

30.3 

1? 

4.958 

9-2 

22.3 

18 

5-283 

9.5 

15.0 

19 

5.607 

5.0 

5.8 

20 

5.930 

9.4 

20.6 

Total  Heat  Transfer  Rate  at  Hot  Plate  for  aT^  =  5-29  Btu/hr 
Total  Heat  Transfer  Rate  at  Hot  Plate  for  AT^  =  9-73  Btu/hr 
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Table  XXI 


Local  Heat  Fluxes  for  2.78 6"  x.  3.00"  Geometry 
Aspect  Ratio  =  0.929,  aT  =  39.?53°F,  *  =  30° 


Position 

Humber 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2) 

1 

0.0 

2.3 

-  2.7 

2 

0.398 

13.5 

-11.8 

3 

0.796 

17.3 

-17.1 

4 

1 . 194 

22.2 

-19.5 

5 

1.592 

20.8 

-18.1 

6 

1.990 

15.4 

-19.0 

7 

2.388 

12.4 

-12.8 

8 

2.  ?86 

4.0 

-  2.0 

Total  Heat  Transfer  Rate  at  Hot  Plate  ~  3 • ^7  Btu/hr 
Total  Heat  Transfer  Rate  at  Cold  Plate  =  -3-33  Btu/hr 
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Table  XX J I 

Local  Heat  Fluxes  for  2.786"  x  3  •00"  Geometry 
Aspect  Ratio  =  O.929,  aT  =  53-569°F,  4-  =  30° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 
Cold  Plate 
(Btu/hr  ft2) 

1 

0.0 

2.9 

-  2.4 

2 

0.398 

41.? 

-  9-1 

3 

0.796 

36.0 

-17.3 

4 

1.194 

25.6 

-21.0 

5 

1-592 

18.4 

-26.2 

6 

1-990 

18.6 

-34.0 

7 

2.388 

11.4 

-43 . 6 

8 

2.786 

3-8 

-  2.0 

Total  Heat  Transfer  Rate  at  Hot  Plate  -  5*4l  Btu/hr 
Total  Heat  Transfer  Rate  at  Cold  Plate  -  -5-09  Btu/hr 


Table  XXIII 


Local  Heat  Fluxes  for  5-9"  x  3 •00"  Geometry 
Aspect  Ratio  =  1.96?,  &T  =  25.035°F,  <t>  =  30° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
( Btu/hr  ft2) 

1 

0.0 

4.5 

2 

0.314 

11.5 

3 

0.628 

16.8 

4 

0.942 

13.1 

5 

1.25? 

14.4 

6 

1-571 

13.2 

7 

1.885 

11.4 

8 

2.199 

7.6 

9 

2.513 

8.4 

10 

2.82? 

6. 6 

11 

3.072 

9.7 

12 

3.386 

10.8 

13 

3-701 

10.8 

14 

4.015 

9-3 

15 

4.329 

7.3 

16 

4.643 

8.6 

17 

4.957 

8.8 

18 

5.271 

5-9 

19 

5-585 

4.7 

20 

5.900 

4.4 

Total  Heat  Transfer  Rate  at  Hot  Plate  =  4.80  Bt 
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Table  XXIV 


Local  Heat  Fluxes  for  5*9"  x  1.5”  Geometry 

Aspect  Ratio=3 -933 .  aT-^13 .435°F,  aT2=39.64°F,  *=30° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aT» 
(Btu/hr  ft2 ) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aT2 
(Btu/hr  f t* ) 

1 

0.0 

2.8 

7-9 

2 

0.318 

12.2 

42.0 

3 

0.635 

13.0 

34.1 

4 

0.952 

8.4 

28.5 

5 

I.270 

7.8 

29.0 

6 

1-587 

9-3 

20.9 

7 

1.905 

7.1 

20.5 

8 

2.222 

.  5.1 

15.8 

9 

2.540 

3.5 

16.2 

10 

2.857 

3-5 

11.5 

11 

3.044 

3.1 

13.0 

12 

3.362 

3.0 

12.0 

13 

3.679 

3.1 

13.5 

14 

3.997 

2.8 

13.9 

15 

4.314 

2.9 

10.5 

16 

4.632 

2.5 

8.5 

17 

4.949 

2.2 

7.7 

18 

5-267 

2.1 

7.0 

19 

5-584 

0.5 

3.3 

20 

5.900 

0.7 

0.6 

Total  Heat  Transfer  Rate  at  Hot  Plate  for  aT-^  =  2.46  Btu/hr 
Total  Heat  Transfer  Rate  at  Hot  Plate  for  aT2  =  8.17  Btu/hr 


Table  XXVI 


Local  Heat  Fluxes  for  5*93"  x  0.525"  Geometry 
Aspect  Ratio=11.295.  aT^W.WF,  at2=3?.869°F,  4-^0° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aT, 
(Btu/hr  ft*) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aTj 
(Btu/hr  ft*) 

1 

0.0 

10.3 

22.6 

2 

0.324 

8.8 

29.6 

3 

0.648 

8.1 

34.3 

4 

0.972 

7.0 

32.6 

5 

1.296 

10.4 

24.5 

6 

1.620 

12.2 

29.7 

7 

1.945 

10.1 

27.2 

8 

2.269 

10.7 

36.8 

9 

2.593 

10.7 

34.8 

10 

2.917 

10.3 

25.1 

11 

3.014 

10.7 

24.8 

12 

3.338 

10.7 

24.9 

13 

3.662 

12.1 

28.6 

14 

3.986 

12.1 

25.3 

15 

4.310 

5-7 

23.7 

16 

4.634 

8.5 

23.1 

17 

4.958 

9-0 

23.1 

18 

5.283 

10.2 

12.4 

19 

5.607 

7.2 

5.7 

20 

5.930 

10.0 

10.6 

Total  Heat  Transfer  Rate  at  Hot  Plate  for  AT-^  =  4.7  Btu/hr 
Total  Heat  Transfer  Rate  at  Hot  Plate  for  AT^  =  12.33  Btu/hr 
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Table  XVII 


Local 

Heat  Fluxes  for  2 

.786"  x  3.00" 

Geometry 

Aspect 

Ratio  =  0.929,  AT  =  25 . 00°F, 

4>  —  0° 

Local  Heat 

Local  Heat 

Distance  Along 

Fluxes  at 

Fluxes  at 

Position 

the  plate  "x" 

Hot  Plate 

Cold  Plate 

Number 

(in) 

(Btu/hr  ft5 ) 

(Btu/hr  ft2) 

1 

0.0 

7.2 

-  2.7 

2 

0.398 

14.1 

-10.5 

3 

O.796 

10.6 

-11.9 

4 

1.194 

6.9 

-10.2 

5 

1.592 

8.7 

-  7.3 

6 

1.990 

11.2 

-11.5 

7 

2.388 

10.4 

-14.7 

8 

2.786 

3.0 

-  6.5 

Total  Heat 

Transfer  Rate  at 

Hot  Plate  --= 

2.22  Btu/hr 

Total  Heat 

Transfer  Rate  at 

Cold  Plate  ~ 

-2.3.5  Btu/hr 

<► 
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Table  XVI II 

Local  Heat  Fluxes  for  2.786"  x  3 •00"  Geometry 
Aspect  Ratio  =  0.929,  AT  =  40.073°F.  *  =  0° 


Local  Heat 

Local  Heat 

Distance  Along 

Fluxes  at 

Fluxes  at 

Position 

the  plate  "x" 

Hot  Plate 

Cold  Plate 

Number 

(in) 

(Btu/hr  ft2) 

(Btu/hr  ft2) 

Total  Heat  Transfer  Rate  at  Hot  Plate  =  3.90  Btu/hr 
Total  Heat  Transfer  Rate  at  Cold  Plate  -  -3.68  Btu/hr 


Table  XXIX 


Local  Heat  Fluxes  for  5-9"  x  3.00”  Geometry 
Aspect  Ratio=l .967,  aT1=22.455°F,  aT2=35- 478°F,  $=0° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  ATi 
(Btu/hr  ft2) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  at2 
(Bxu/hr  ft2 ) 

1 

0.0 

3-4 

2.4 

2 

0.314 

10.3 

15.1 

3 

0.628 

10.6 

11-5 

4 

0.9^2 

6.0 

13.2 

5 

1.257 

9.1 

8.0 

6 

1.571 

9.1 

4.0 

7 

1.885 

10.1 

4.0 

8 

2.199 

9.2 

4.0 

9 

2.513 

9.2 

8.0 

10 

2.827 

6.0 

13.2 

11 

3.072 

?-5 

19.0 

12 

3.386 

8.3 

20.6 

13 

3.701 

9.4 

20.4 

14 

4.015 

11.7 

22.5 

15 

4.329 

11.8 

22.0 

16 

4.643 

9.4 

20.7 

17 

4.957 

4.6 

19.2 

18 

5.271 

8.8 

10.5 

19 

-  5.585 

8.2 

8.2 

20 

5.900 

4.4 

11.3 

Total  Heat 

Transfer  Rate  at 

Hot  Plate  for  AT-^ 

=4.28  Btu/hr 

Total  Heat 

Transfer  Rate  at 

Hot  Plate  for  ATg 

=  6.55  Btu/hr 

Table  XXX 


Local  Heat  Fluxes  for  5*9"  x  1.5"  Geometry 
Aspect  Ratio  =  3-933,  aT  =  21.8l8°F,  *  =  0° 


Position 

Number 

Distance  Along 
the  plate  "x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
(Btu/hr  f t2 ) 

1 

0.0 

4.9 

2 

0.318 

16.3 

3 

0.635 

13-9 

4 

0.952 

8.1 

5 

1.270 

9-0 

6 

1.587 

8.2 

7 

1.905 

6.2 

8 

2.222 

6.0 

9 

2.540 

6.0 

10 

2.857 

6.1 

11 

3.044 

7.4 

12 

3.362 

7.4 

13 

3.679 

7-7 

14 

3.997 

9-1 

15 

4.314 

10.3 

16 

4.632 

11.1 

17 

4.949 

10.0 

18 

5.267 

7.4 

19 

5-584 

6.3 

20 

5.900 

4.2 

Total  Heat  Transfer  Rate  at  Hot  Plate  =  4.21  Btu/hr 
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Table  XXXI 


Local  Heat  Fluxes  for  5-9"  x  0.75"  Geometry 
Aspect  Ratio=7. 86?,  aT-l=18.803oF,  aT2=38.304°F,  *=0° 


Position 

Number 

Distance  Along 
the  plate  ”x" 
(in) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  AT! 
(Btu/hr  ft2 ) 

Local  Heat 
Fluxes  at 

Hot  Plate 
for  aTj 
(Btu/hr  ft2) 

1 

0.0 

1.9 

2.4 

2 

0.315 

12.7 

22.0 

3 

0.630 

9.6 

26.? 

4 

0.944 

10.1 

24.2 

5 

1.259 

11.7 

27.0 

6 

1.574 

9-3 

19.1 

7 

1.889 

9.0 

18.4 

8 

2.204 

9-9 

18. 9 

9 

2.518 

9.3 

24.9 

10 

2.833 

9.1 

19.6 

11 

3.066 

9-5 

20.0 

12 

3-381 

10.2 

12.8 

13 

3-696 

9.8 

18.5 

14 

4.011 

6.6 

18.9 

15 

4.325 

7.9 

18.1 

16 

4.640 

8.7 

19.5 

17 

4.955 

9.1 

20.7 

18 

5.270 

8.3 

22.2 

19 

5-585 

5.2 

13.5 

20 

5.900 

4.0 

5.4 

Total  Heat  Transfer  Rate  at  Hot  Plate  for  AT^  =  4.36  Btu/hr 
Total  Heat  Transfer  Rate  at  Hot  Plate  for  aT2  =  9-52  Btu/hr 


Table  XXX1J 


Local  Heat  Fluxes  for  2.825"  x  3*05"  Geometry 
With  Diagonal  Partition 
Aspect  Ratio  =  0.926,  at  =  39-355°F,  *  =  90° 


Pcsi 

Num 

tion 

ber 

Distance  Along 
the  plate  "x" 

(in) 

Local  Heat 
Fluxes  at 
Hot  Plate 
( Btu/hr  ft’) 

Local  Heat 

Fluxes  at 

Cold  Plate 
(Btu/hr  ft2) 

1 

0.0 

0.9 

0.0 

2 

0.327 

18.8 

-  4.0 

3 

0.653 

15.8 

-  2.4 

4 

0.980 

10.6 

-  5*6 

5 

1.306 

9.6 

-  9-2 

6 

1.519 

8.6 

-10.5 

7 

1.845 

7.1 

-12.1 

8 

2.172 

1.8 

-13.9 

9 

2.498 

6.1 

-18.0 

1 

0 

2.825 

0.0 

-  0.9 

Tote 

il  Heat 

Transfer 

Rate  at 

Hot  Plate  - 

2.06  Btu/hr 

Tot; 

ll  Heat 

Transfer 

Rate  at 

Cold  Plate  = 

-I.98  Btu/hr 

1 
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setting.  The  average  Nusselt  number  was  found  to  be  primarily  a 
function  of  Grashof  number  or  Rayleigh  number.  Three  correlation 
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equations  developed  by  earlier  studies  are  fitted  to  the  data. 
The  correlation  equations  are  found  to  approximate  the  heat  flux 
within  +  11$,  up  to  a  Grashof  number  of  1.2  million. 
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